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Abstract 

The dynamical response of a lipid membrane to a local perturbation of its molecular symmetry 
is investigated theoretically. A density asymmetry between the two membrane leaflets is predom- 
inantly released by in-plane lipid diffusion or membrane curvature, depending upon the spatial 
extent of the perturbation. It may result in the formation of non-equilibrium structures (buds), 
for which a dynamical size selection is observed. A preferred size in the fim range is predicted, 
as a signature of the crossover between membrane and solvent dominated dynamical membrane 
response. 
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The formation of small membrane structures (vesicles or tubules...) is required for most 
inter and intra cellular transports in biological cells Q]. While important progress have been 
made in the identification of key membrane proteins, recent work have focused on the lipid 
molecules themselves Q]. Of particular interest is the lipid translocation by specific enzymes 
(flippases), and the membrane morphological changes they trigger P, ^\. From a physical 
point of view, the formation of a bud from a fluid membrane has mostly been considered 
as the result of phase separation in mixed membranes j^], or as a global shape transition 
of closed membranes due to geometric frustration The final membrane conformation 
then corresponds to a global equilibrium, and is not expected to strongly depend upon 
dynamics (7]. Small invaginations of the plasma membrane of biological cells should however 
result from localized perturbations, rather than from global changes at the scale of the cell. 
This paper present a theory for the dynamics of relaxation of one such local perturbation; 
a local membrane asymmetry, created by a sudden flip of a number of lipid from one leaflet 
of the bilayer to the other. This situation is of fundamental interest, as it may result in 
the formation of non-equilibrium membrane structures. The model also seeks to mimic 
experiments where the local perturbation of a giant vesicle (~ 100/im) by adsorption of 
DNA molecules results in the production of small dynamical vesicles (~ /im) with their 
membrane packed with adsorbed DNA [8]. 

In the present theoretical analysis, it is shown that the local membrane perturbation 
introduced by lipid translocation (stretching of the depleted monolayer, and compression of 
the enriched one, Fig[TJ), may be released both by opposite monolayer flows, leading to a 
diffusion of the perturbation (FigOJi) or by membrane curvature (FigHt). While the former 
mechanism potentially leads to a vanishingly small energy, the faster of the two processes will 
control the membrane relaxation. It is shown below that large-scale membrane curvature 
is hindered by bulk flow, while small scale-membrane curvature is slowed by membrane 
flow. The cross-over between these two regimes defines a critical dynamical lengthscale 
Xd ~ 1/ini in water and ~ lOOnm in the more viscous inner cellular environment, at which 
the perturbation is optimally converted into transient membrane curvature, and fully formed 
membrane bud are most likely to be observed. 

Lipid membranes are self-assembled fluid bilayers. Their equilibrium properties such as 
lipid density O and membrane thickness 2d (~ 4nm), result from a balance of the hydropho- 
bic attractions between the lipid tails by steric or . Deviation from 
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the optimal density costs an elastic energy (per monolayer, per unit area) K s ((p — o ) 2 /(20o)? 
with a large stretching modulus K s ~ 2hksT '/nm 2 10( (ksT is the thermal energy). The 
membrane bending energy kC 2 /2 involves the local membrane curvature C and a fairly small 
bending modulus k ~ 25ksT. Noting the outer (+) and inner (-) monolayer densities (at 
the mid-plane) ij)± (with <f>± = ip±(l =F dC)), the membrane elastic energy is best expressed 
in terms of the local average dilation p = t±+?t= _ \ anc i the dilation difference between the 

: 

Felast = J dS(K s (p - dCf + K s p 2 + ^C 2 ) (1) 

Bending a symmetrical bilayer (p = 0) involves a bending modulus k + 2K s d 2 ~ 200A;bT. It 
is reduced to k in fluid membranes since the monolayer densities can adjust to p = dC. The 
(small) "bending ratio" of these two quantities e = K,/(2K s d 2 ) will prove important later. 

Here, we study the membrane response to the sudden establishment of a local concen- 
tration asymmetry: some fraction po of lipid being flipped from the "down" to the "up" 
monolayer over an area Sq (Fig^). For simplicity the membrane bears no global surface 
tension, which is to say that p = throughout. Lipid flow toward the stretched region of 
the "down" monolayer and away from the compressed region of the "up" monolayer, diffuse 
the perturbation over the membrane (Fig^,). Concomitantly, the compression-extension 
creates a bending torque leading to membrane curvature (FigO^)- Diffusion of the pertur- 
bation is limited by the sliding friction between monolayers, while membrane curvature is 
limited by membrane and solvent viscous dissipation. 

Local balance between elastic and dissipative forces yields local equations for the evolution 
of the membrane shape. Such equations exist for small membrane deformations|2|, but are 
hardly tractable for large deformations. In consequence, we take a simpler approach where 
the perturbed membrane (of area S(t)) is parameterized by a spherical cap of constant 
curvature C(t) and constant dilation difference p{t), connected to an unperturbed (flat - 
C = p = 0) membrane, see Fig^. This approximation disregards the existence of a 
connecting neck between the bud and the flat membrane, the stability of which probably 
influences the scission of the bud by membrane fusion |l2j]. In this work, buds are simply 
assumed to detach beyond a critical value of the budding parameter (defined in Eq.((2|)). 
In the biological context, it may correspond to a neck small enough for membrane proteins 



such as dynamin to pinch off the bud 



3 



Spontaneous lipid flip-flop between the two monolayers is very unlikely over the evolution 
time, which imposes the conservation of asymmetry S(t)p(t) = SqPo- The stretching stress 
in the perturbed region is released when C = p/d (Eq.([T|)). The membrane geometry is 
described by a budding parameter, equal to unity for closed buds: 

< B " s w < B ™ (2) 

where B d °' = Sopl/^Qird 2 ) corresponds to a situation where the perturbation is entirely 
converted into curvature, without diffusion along the membrane. As we will see, Bd is 
always much smaller than B d °\ The optimal conversion of perturbation into curvature is 
observed for a particular dynamical lengthscale \/So ~ Ad (Eq-©)- 

Dynamical equations for the two independent variables S and C are obtained using a 



Lagrangian description 14], where elastic and dissipative "forces" are calculated from the 
variation of the elastic energy J r e iast[{P: C}}, and the energy dissipated per unit time Vdiss'- 

d{p^c} + d{d tP ,d t c} - [ ' 

Variation with p accounts for the migration of lipid molecules under a gradient of chemical 
potential, and variation with C describes the membrane deformation due to the bending 
torque. 

The three sources of dissipation involve three constitutive parameters. The viscosi- 
ties of the solvent 77 (~ 10~ 3 N/m 2 for water) and of the membrane p m (~ 10 -9 N.s/m), 
couple to gradients of bulk and membrane velocity fields v into viscous shear stresses 
crjj = 77 (djVi + diVj) [lq . An interlayer friction parameter b m (~ 10 8 N.s/m 3 ), couples to 
the velocity difference 5v between membrane the two monolayers in relative motion[10j. 
Note that p m and b m are related to the membrane (3-dimensional) viscosity r] m by the scal- 
ing relations r] m ~ p m /d ~ b m d, where the membrane viscosity is typically a thousand times 
the viscosity of water r) m ~ 1 N/m 2 [l^|. The total dissipation is Vdiss = Vb m + + V v , 
with 



TV = ^ / dS5v 2 ■ V 



l -jdV^f (4) 

ij 



2 J 7 ' " 2 V 

The membrane viscous dissipation V^ m is obtained by substituting p m to r\ in V v and 
integrating over the membrane surface instead of the volume. 

The three contributions to the energy dissipation are calculated as follows (see for 
a related calculation). The bud volume V and neck aperture L satisfy V = (l — §-Bcz), 



ttL 2 = S(l — Bd). The curvilinear, radial, and spherical coordinate systems defined FigO 
are used to parameterize flows in the curved and flat part of the membrane, and in the 
bulk, respectively. The differential monolayer velocity Sv is present in the perturbed part 
of the membrane. It is axisymmetrical, and relaxes the concentration difference p according 
to the continuity relation^ V c .<5v = —d t p, or 5v = — 4p(l — cos ip) / (C sin ijj) (where the 
dot represents a time derivative, V c is the curvilinear gradient along the membrane, and 
ip is defined Fig|2^. The intermonolayer friction part of energy dissipated per unit time is 
calculated from Eq.(JU): 

^sM^K^T*)) ( 5) 

We assume incompressible flows in the membrane and the surrounding fluid. Most of the 
membrane flow comes from bringing membrane area from the flat membrane into the bud. 
The flow matches the variation of the area AS = S — nL 2 = SB^, and corresponds to a 
radial velocity field tWmbmne = — dt(AS) / (2nl) (where / is defined FigEb)- The main source 
of solvent flow comes from the variation V of the volume inside the bud, which imposes a 
flow going through the circular neck of radius L. In the spherical coordinate system FigOJ, 
the velocity field is radial and have the expression [3]: = ^3 cos 2 6. The bulk and 

membrane viscous dissipation are obtained from Eq. (|^jl 

Viscous dissipation in the solvent occurs in both sides of the membrane. V v in Eq.© is 
valid for the fluid entering the bud. For the fluid expelled by the bud, the neck size L is 
approximately replaced by the bud radius 2/C beyond hemispherical bud. Extra sources of 
dissipation, such as the solvent flow inside the bud, and membrane shear flow in the curving 
membrane, have been checked not to modify the membrane relaxation qualitatively, and 
are omitted here for clarity. The evolution of S and C is entirely determined by Eqs(Q3 01 
EJ El)- Notwithstanding the simple assumptions for the membrane shape, these equations 
(Eqs. (jl()pil)) ) are fairly complex and have to be dealt with numerically. In the limit of small 
deformations Bd 1, all equations may be linearized, and membrane viscous dissipation is 
negligible. The membrane asymmetry relaxes via diffusion along the membrane according 
to the expression Sp = —D(p — dC) with the diffusion coefficient D and characteristic 
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relaxation time U 



0.1ms for So ~ fim 2 ): 
K 



D = 4:71- 



t n = 



15 



(7) 



The influence of the perturbation lengthscale Sq is emphasized by using dimensionless vari- 
ables for area, curvature and time: 



S_ 

So 



dC_ 



t 

T= t~Q 



(8) 



The new variables are bounded by: 1 < s < oo and < c < 1, and the budding parameter 
is Bd = B d °^ sc 2 (see Eq.(j2J). Balancing the curvature forces lead to an equation for C. 



The interplay 
X D ~ bd 2 /rj§, 



aetween membrane and solvent dissipation defines a characteristic lengthscale 



y 



. to which the size of the perturbation should be compared. 



2O0F bd 2 



S n = 



Sq 



(9) 



11 V Xd 

This lengthscale A^> is of order 1/xm in water, but is much shorter (~ lOOnm) in biological 
condition, since the cytosol can be very viscous. 

Expressed in the dimensionless variables, the set of equations reads[18]: 



1 



B, 



1 



sc) 



(10) 



/7r(l-25 d ) 2 / (l-B d ) , /9 V 

+ fe)(r^k)( s2 ' i+ ^)=( 1 -( 1+ ^ <"> 

The driving forces for diffusion and membrane deformation (RHS of the above equations) 
both show the competition between increases of s and c. 

At short time (s ~ 1 and c <C 1), the evolution is given by s ~ 1 and c ~ 1/a/^o- 
For So < 1, membrane curvature moves little solvent and occurs faster than diffusion. 
On the other hand, membrane deformation is slow for Sq > 1. For larger deformation 
however, the membrane viscosity comes into play. It is characterized by a single parameter 
a m = fi m /{b m d 2 ) (typically slightly larger than unitvfl^|). and becomes important for small 
scale perturbations (Bd > ySo/a m )- 

Eqs. (jlOHH) show that membrane deformation will evolve non-monotonously. Indeed, the 
driving force for diffusion ((1 — sc) - Eq. fjlOj) ) vanishes for complete release of the stretching 
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stress (sc = 1, or C = p/d). On the other hand, the driving force for membrane curvature 
((1 — (1 + e)sc) - Eq. fllljO vanishes for smaller deformation because of membrane bending 
rigidity (e = n/(2K s d 2 )). As a consequence, the budding parameter rises from zero to 
a maximum value B d max ^ over a time t grmnt h, then it decays slowly to zero over a larger 
time tdecay The value of the maximum is a measure of the amount of initial perturbation 
converted into curvature. It will be discussed extensively in the remaining of this letter, for 
it determines whether the blister may turn into a well-formed bud before the perturbation 
diffuses out. For small deformations, the rise occurs over the linear evolution time t growth ~ 
to, and the decay, mostly driven by bending rigidity, is of order tdecay ~ to/e. These times are 
however much larger in the non-linear regime. To give a feel for the numbers; for So = 1, the 
maximum deformation is an hemispherical buds (B d max ^ = 0.5) for an initial perturbation 
Bjp = 3, while nearly completed buds ^B d max ^ = q.9) require B d > 6 (see Fig|3J). Fixing 
Xd = lpm, hemispherical buds have evolution times t growt h = 0.7ms and tdecay = 26ms, and 
nearly completed buds t growt h = 7ms and tdecay = 300ms. 

Many results can be extracted from Eqs. ()10Jllj) . which determines (within the approxima- 
tions) the full dynamics of the membrane deformation. Here we concentrate on the following 
question. How much should we perturb the membrane in order to see buds? FigJS] shows 
the value of the initial density of flips, and of the initial perturbation B d °^ needed to obtain 
a given maximum of the budding parameter B d , as a function of Sq. One may practically 
picture B d as some critical value at which the neck connecting the bud to the mother mem- 
brane becomes unstable. More generally, it illustrates the non-trivial competition between 
membrane and solvent dynamics for large membrane deformation. 

FiglHl shows that membrane viscosity strongly reduces the likeliness of large membrane 
deformations (B d > 0.5) for small initial extensions So < 1. More extended perturbations on 
the other hand, are dominated by solvent dissipation. Interestingly in this case, the diffusion 
of the perturbation (increase of S) promotes bud closure (increase of C), to minimize the 
variation of the bud volume (V = when SC{2Bd — 1) = 2CS{\ — Bd))- In consequence, a 
bud is likely to close if it reaches the hemispherical shape, and the results in FigO overlap 
for Sq ^> 1. Furthermore, the density needed to reach a given membrane deformation B* d 
surprisingly tends to a constant value, so the total perturbation mass p So linearly increase, 
as the extension S increases. 

The interplay between membrane and solvent dissipations produces a minimum of B d 
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for a given Sq of order A^ 2 , purely dictated by dynamics. Mature buds have best chances to 
be produced at this particular lengthscale, as the maximum membrane deformation B* d is 
closest to the optimal one B d °\ Interestingly, the DNA-induced budding of giant vesicles 
results in endosomes of size close to this optimal dynamical size in water (~ /im). It is 
hazardous to make more quantitative comparison with these experiments, since a number 
of other factors such as the kinetics of DNA adsorption certainly influences the formation 
of buds, but it is one prediction of the present work that the bud size can be altered by 
changing the solvent viscosity. 

According to FigEl well- formed buds with Bd = 0.9 require flipping only p ~ 3.5% lipids 
over an area Sq ~ 1/zm 2 (for = 1/im), while it necessitates po — 35% for Sq ~ (200nm) 2 . 
This amounts to flipping about 10 5 molecules in both cases. The flips need not being 
perfectly synchronized, as successive events may act cooperatively if they are separated by 
a time lag shorter than the perturbation decay time t decay As already discussed, t decay ~ 
to/e ~ So/(De) in the limit of small membrane deformation. Note that the evolution becomes 
much slower in the large deformation regime, where dissipation is very large. 

In conclusion, non-linear dynamical equations for the formation of bud-like invaginations 
in fluid membranes have been derived, and applied to the relaxation of a localized pertur- 
bation of the up/down symmetry of lipid bilayers. Non-equilibrium buds may be produced 
provided that relaxation by membrane curvature occurs faster than the diffusion of the 
perturbation along the membrane. The membrane response critically depends upon the 
extension of the initial perturbation, which determines whether the membrane deformation 
is limited by membrane or solvent dynamics. Bud formation is most likely at the crossover 
between the two regimes. 

I thank Albert Johner for many enjoyable discussions. 
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FIG. 1: Relaxation of a membrane asymmetry (a) in-plane diffusion and (b) membrane curvature. 
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FIG. 2: Parameterization for the membrane shape surface S, curvature C = 2/R and radius of the 
neck L, and the various coordinate system used to derive the dissipation 
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FIG. 3: Density of flips po and strength of the initial perturbation B d needed to reach hemi- 
spherical and nearly-completed buds, as a function of the initial extension So (e = 0.2 - = 3 - 
Xd = 1/um). Dashed lines show the densities po corresponding to = 0.5 (grey) and 0.9 (black). 
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